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Abstract 

Mellin transform is used to evaluate an integral involving the product of four Bessel functions 
and a power. Using this method the result is obtained in terms of generalized hypergeometric 
functions 6-^5- 

PACS numbers: 33C10 



I. INTRODUCTION 



The integral that we want to solve in this paper is of the type: 

POO 

/ dx x^J a (ax)Jfs(ax)J^(bx)Js(bx). (1) 
Jo 

In our analysis the two cases: a > b and b > a will be discussed. Our intention is 
to evaluate the above integral using the Mellin transform for representing relation ([I]) as 
a contour integral in the complex plane, whose integrand will contain ratios of gamma 
functions. The evaluation that will be presented here follow the method used in based 
on Mellin transforms, that can be used for evaluating the integrals of this type. 

The paper is organized as follows: in Section 2 we introduce the Mellin transform for a 
Bessel function as in [2] and we write the integral as a contour integral in the complex plane. 
In Section 3 we present the main steps that help us to arrive at the final result. In Section 
4 we point out the possible applications in the scattering theory on an expanding universe. 



II. MELLIN TRANSFORM 



The integral ([T]) can be brought in the form 



o 



where we use the notations rj = ax and r = b/a. 

The solution that we obtain is valid for a range of parameters consistent with the conver- 
gence of the integral from relation (J2]). A study of this expression reveals that the integral 
converges when: 

Re(n) < 2, 

Re(fi + a + (3 + 7 + 5 + 1) > (3) 



where the first expression is obtained by considering the behavior of the integrand of (j2J) as 
rj approaches infinity and the second expression is obtained by considering the behavior of 
the integrand as rj approaches to zero 

The integral that interest us can be brought in the following form: 



oo 



( /U ,r)= / dr ] ^M 1 (r ] )M 2 (r V ) (4) 



where M-y{rj) and M2{rrj) are functions whose Mellin transform can be obtained as ratios 
of products of gamma functions. Let us start by considering the Mellin transform of (j2J). 
The Mellin transform of a function M{rj) will be in our notations A4(s), and will have the 
following expression: 

POO 

M(s) = / dT]Ti s - x M{rj). (5) 
Jo 

With the observation that the original function is recovered from its Mellin transform eval- 
uating the contour integral: 

i pX+ioo 

M( V ) = — ds V - s M(s), (6) 

Zm JX-ioo 

where A is a small positive constant. The Mellin transform of Ia^si^, t) will be [2[: 

poo 

T aj3lS (ii,s) = / dTT s ~ l I a(3l& (n,T). (7) 
Jo 

Now substitute (BJ in equation (C3) and using equation ([5]) we obtain: 

Ia/3 7 s(n,s) =M 1 (n-s + l)M 2 (s). (8) 

It not difficult now to use the Mellin inversion formula given by (EI), for obtain an alternative 
expression for (T5]): 

■y rX+ioo 

Iap 1 s{n,r) = — I dsT~ s Mi{{i- s + l)M 2 {s). (9) 

2m JX-LOO 

In our case the functions whose Mellin transform can be expressed as ratios of gamma 
functions are: 

Mi = J a (v)Mv), 

M 2 = J^JsiTT]). (10) 

The required Mellin transforms are: 
Mi(n - s + 1) - 



r(i - g)r( 7+ ^ +s ) 
-^ 2 ( s ) ~~ 2 1-s r ( ~^ +<5 ~ s + 2 ) y ( t +,5 - s + 2 ) r ( ~<- s - s+2 ) ' 

Replacing (ITTj) in we obtain: 



1 pA+tOO 



2 ,s ~'"r( ~° M )r( M )r 
r(i- g) r(^) 2 2 



Making the coordinate transformations s = 2u and c = A/2 as in 2] and using the 

n 

duplication formula for gamma functions [1] 

T(2z) = (2ir)- l/2 2 2z - 1/2 T(z)T(z + 1/2) (13) 

we obtain the following result: 

—n f c+io ° T(v - ^)T(v + i±±l)r(^±M±E±l _ v ) 
W^r) = -1 / dur- 2u 1^ 2 )t{u+ 2 M 2 «j 

w ; 47r 2 y c _ loo r( -°+w +u)r( a+/? - M+ i +M )r( °-w +lt ) 

r(i- M )r(i- M )r(^ + M ) 
r(^±2 _ M)r( 7±|±2_ M)r( 7z|±2 _ u y ^ > 

The above integral is a contour integral in complex plane. As one know the poles of the 
gamma function occur when their arguments are zero or negative integers. Then the integral 
can be evaluated using the fact that the poles of the numerator appears for the above 
mentioned values. Since there are an infinite number of such poles inside the contour an 
infinite series will be obtained when the residue of these poles are summed. 



III. EVALUATION OF THE INTEGRAL 

The integral from relation fl!4p will be evaluated using the complex-variable technique. 
Each gamma function from numerator has a pole when the argument equals zero or is an 
negative integer. The poles of the integrand are: 

u = — — k: u = — k ; 

2 2 

„__*_(!+« (15) 

and 

u = 1 + k; u = — h k ; 
fi + a + fi + l 

u= + k, (16) 

with k = 0,1, 2, 3.... For the following calculations we assume that the values of parameters 
H, a, j3, 7, 5 are such that a value for c can be found for which the poles from expression (fl5l) 
lie to the left of line u = c while the poles given by ( fTBT) lie to the right of line u = c. 
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Let us begin with the case a > b (r < 1), when the poles represented by (IT51) give 
contributions. For evaluating the integral (|14p we use the following relation: 



00 (-l) k 

dzH{z)T{z) = 2m V ^—^H(-k), (17) 
c 'i k=o kl 

where H(z) is a analytic function without singularities in the left half-plane. Using ( ITT)) we 
find that the final result is: 



L S fc! r( z^±i _ A;)r(^±f±i - ifc)r(2=p - fc) 



X 



r( -^-f+ 2 + fc)r(^^ + fc)r(^^ + k) + ^ fc! - fcj 

r(-| - fc)r(fc + ^l)T(9±m + fcjrgg + AQIX 7 ^- 1 - fc) 

x r(^ - fc)r(s=£ - fc )r( -7+^+3 + ^ 7+^+3 + fc ^ r ^ 7 -^+ 3 + ^ 

~ (_1)* r 7+5+2/= r (7±|±i + jfe) r (21a±£ + jfe) r (z2=5zfi _Jfe) 

+ ^ fc! r ( °+/3-M-7-«?+l _ _ fc)r(^±£ -M-7-«?+l _ fe ) r ( a-/3- M -7-^+ l _ fc ) 



X- 



fc=0 V 2 'V-*- V 2 ' J - 2 

-7+6-/^+1 _ ^■jp^ 7+»5+a+/3+M+l _|_ ^ 



(18) 



T( 7 + + i)r(5 + fc + i)r( 7 + 5 + k + 1). 

This result can be written now in terms of generalized hypergeometric functions, using the 
following relation: 

f, z k n£ r(fc + go ns r(m t - fc) n£i iw nS 

S fc! n£i r(fc + n£ rfe - fc) n£i r( ft ) n2i rfo) 

N g +N P F Nm+Nq {g u g 2 , ...,gN g , -pi + 1, -P2 + 1, ••• -Ptv p + 1; 91,92, 

9^, -m a + 1, -m 2 + 2, -m^ + 1; (-1)^"^). (19) 

Using equation (fT9l and the following relations 

r(i/2) = tt 1 / 2 , r(-i/2) = -27T 1 / 2 , 

r(*)r(z + 1/2) = 2 1 - 2 % 1/2 r(2^) (20) 
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the final form of equation ([TBI) is: 

y-l+H T-y. Y(— 



2-i+n T -» r(^ti)r(^±|±i)r(2±|±M) 

Q+ff+l \ri/a- 3+1 \r/ \p/ ~7+<S— /j+2 ^Tn/7+<5-/j+2 ^ 



1 /a + /? + 1 l- n 2-fi a- (3 + 1 -a- (3 + 1 

X r ^-^-M+2 y 6 5 ^ 2 ' 2 ' 2 ' 2 ' 2 : 

-a + (3 + 1 -7 + 5 - /a + 2 7 + 5- /^ + 2 7-5-/1 + 2 1 



_ 7 _ 5 _ /i + 2 _ 2 \ 2- 2 +^r-^+ 1 r(2-/i)r(^f^)r( 7+ ^ 1 - 



2 

+ 

2 



r(^)r(^)r(^)r( - 7+ 2 " At+3 )r( 7+ ' 5 7 +3 ; 

a + (3 + 2 2-/i 3-/i a-/3 + 2 -a - /? + 2 



X p^ 7-(5- M +3 ^ 6^5 

-a + /3 + 2_-7 + 5 - /a + 3 7 + 5 - /a + 3 7-5-/1 + 3 3 
2 ' 2 ' 2 ' 2 ' 2' 

_ 7 _ 5-/X + 3 2 \ 2^r^r(^^)r( Q+ ^ 2 +7+m )r(^±|±^) 
^ ;t H " 



p/ a+/3— 7— 5— jj+1 a— /3— 7— <5— /x+1 -a+/3— 7— <5— /j+1 ■ 

1 /a + /5 + /i + 7 + 5 + l 7 + 5 + 1 

6-^5 ' 



r( 7 + i)r(5 + i)r( 7 + 5 + i) u J V 2 

7 + 5 + 2 a-/3 + /i + 7 + 5 + l -a-/3 + /i + 7 + 5 + l 

2 2 ' 2 

-« + /? + // + 7 + 5+1 r 1 -I , r , 1 I 1 /i + 7 + ^ + 2 

;<J +1,7 + 5 + 1,7 + 1, . 



(21) 



2 

/i + 7 + 5 + 1 



When r > 1 the contour is completed in the right half-pane. In this case the poles that 
give contributions are given in equation (1T61) . Now the sign of equation (I14p must be changed 
because the contour is traversed in the opposite direction. The contribution due to the poles 
will be evaluated this time using: 

\k 



r (— i) k 

j> dzH{z)V(-z) = -2iri22^-H(k), (22) 
where H(z) is a analytic function free of singularities in the right half-plane. Using (l2~2~l) to 



6 



evaluate tHM we finally obtain: 



I a /3y6(fJ>,T) = — 



■ ~ (_!)* r -(«+,3+M+2fc+l) + 0+1+1^(0+1+2 + q 

L ^ fc! r(/? + l + fc)r(a + /? + l + fc)r(a + l + fc) 



r(^+^ - fc)r( 1 -"/^ - fc)r( 7+,5+Q + /3+;j+1 + k) ^(-i)' 
x r( -t+^-q-^-m+i _ /;)r( 7+<5 ~"~ /3 ~ A ' +1 — /c)r( 7 ~ <5 ~"~ /3 ~ At+1 — &;) fc! 



2 'V- V 2 'V- V 2 'V fe=0 

r -2fc-i r(i _M + fc)r (i _ fc )r(«+|±tt _ k)v(^ + fc) 



r( ^-" +/3+2 + fc)r( a+/3 7+ 2 + fc)r( "- /3 7 +2 +)r(^±i - fc)rp±|±i - *) 



2 1 V 2 1 V 2 

X 



X 



r(2±|±i + k) ^ (-i) fe r~ 2 - 2fc r( Q+ y- 1 - fc)r( 2 ^ + fe)r(^ + fc) 
r(2z|±i _ fc) + fci r( a+/3 7+ 3 + fc)r( ~ Q+/3 2 " M+3 +~fcj 

r(^ + fc)r(-fc-|) " 

r( Q - /3 ~ M+3 + k)v{=^- - fc)r(2M - fc)r(2=^ - fc). ' 



2 1 '"J 1 - \ 2 '"/^ V 2 "V x V 2 

For obtaining the final result we must use now (120]) and (|T9|) : 

2/i t -(q+/3+m+1) p( -Q-^~M ^p( Q+/3+1 \p/ a+/3+7+(5+^+l ■ 

i> + i)r(/? + i)i> + p + i)r( -^-"-^+ I)f^ ; 



Ia/3-ys(f J ', T ) - > -i -.Tn/ -7+5-Q-/3-M+l -|-p( 7+^~a-/'-;'-H 

'« + /? + ! a + /? + 2 7 + <5 + a + /? + /i + l 



~2 



-p( 7-5-a-/3-M+l ' 

7~5 + a + /? + /i + l -7-5 + a + /3 + /i + l 

2 ' 2 ' 

-7 + 5 + a + /? + /i + l « + /? + /i + 2 
; p + 1, 5 + p + 1, a + 1, 



2 

a + /3 + // + i _ 2 \ 2 - 1 +^r- 1 r(i^)r(^^)r(2±|±i; 



r 



2 ' y r( a+ ^ +2 )r( Q -^ +2 )r(^^)r(^^)r(^|±^ 

1 (l-n 2-fj, 7 + 5 + 1 7-5+1 -7-5+1 

X r( /3 ~ a ~ At+2 ) 5 ' ~~ 2~~ ' 2 ' 2 ' 2 ' 

-7 + 5 + 1 -a + /3-/i + 2 a + /3-/i + 2 a - /3 - /i + 2 



-a - (3 - /i + 2 i _ 2 \ 2 - 2 +^r- 2 r( 2 ^)r( a+/3 ^- 1 )r( , 



X 



T+tf++ 

_. -2 \ ' V~; x V 2 ^ V 

2 '2' / p( Q+/3-M+3 -jp( a-/3-^t+3 -jp( -Q+/3-/J+3 \ 

2-/i 3-/i 7 + 5 + 2 7-5 + 2 2 - 7 -5 



r(izj)r(^)r(~ 7 2 + ' 5 ) V 2 ' 2 ' 2 ' 2 ' 2 

7 + 5 + 2 /3 — a — ji + 3 a + P — fi + 3 a — [3 — fi + 3 —a — (3 — \i + 3 



2 

2' 



This is the result of the integral in the case r > 1. 



IV. CONCLUSION 



We presented here a method based on Mellin transform to evolve one integral that contain 
the product of four Bessel functions with two different arguments. The integrals presented 
here could be of interest in the perturbative QED on de Sitter space and in the study of 
Feynman propagators. It is known that in de Sitter case the solutions of the free field 
equations for Dirac and scalar fields are given in terms of Hankel functions which can be 
expressed with the help of Bessel function. Also in the present literature the scattering 
amplitudes on curved space seems to receive a little attention. An explicit calculation show 
that the scattering amplitudes in the first order of perturbation theory, imply integrals with 
products of two Bessel functions. Hoverer in second and fourth order of the perturbation 
theory the situation seems to be more complicated with integrals that contain products of 
multiple Bessel functions. 
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